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Fragmentation and quark distribution functions for the pion and kaon 
with exphcit flavor- SU(3)-symmetry breaking 



We investigate the unpolarized pion and kaon fragmentation functions, employing the nonlo- 
cal chiral-quark model, which manifests the nonlocal interactions between the quarks and pseu- 
doscalar mesons, considering the explicit flavor-SU(3)-symmetry breaking in terms of the current- 
quark masses. Moreover, we study the quark-distribution functions, derived from the fragmentation 
ones with the Drell-Yan-Levi relation. Numerical results are evaluated to higher by the DGLAP 
evolution and compared with the empirical data. The ratios between the relevant valance quark- 
distribution functions are also discussed. It turns out that the present results are in relatively good 
agreement with available data and other theoretical estimations. 

PACS numbers: 12.38.Lg, 13.87.Fh, 12.39.Fe, 14.40.-n, ll.lO.Hi. 

Keywords: Kaon and pion fragmentation and quark-distribution functions, flavor- SU (3) -symmetry breaking, 
nonlocal chiral-quark model, DGLAP evolution. 



To apply perturbative Quantum Chromodynamics (QCD) to study hadronic processes such as the deep-inelastic 
electron-proton scattering or semi-inclusive hadron production process, provided that QCD factorization theorem is 
applicable there, the cross section is able to be expressed as a convolution of the two parts: the process-dependent 
perturbative QCD(pQCD) calculable short-distance parton cross section, and the universal long-distance functions 
which can be extracted from experiments, but cannot be calculated by pQCD because the strong interaction in the 
long distance is nonperturbative. 

These long-distance functions are the fundamental nonperturbative ingredients for the analysis of the scattering 
processes involving hadrons. The fragmentation functions and the quark-distribution functions are belonged to this 
category. The fragmentation functions, Dq{z) indicate the probability for a hadron fragmented from a quark with 
the momentum fraction where the subscripts h and q indicate the hadron and quark, respectively. It plays 
an important role in analyzing the semi-inclusive processes in the electron-positron scattering, deep-inelastic proton- 
proton scattering, and so on [1 11 . The parton distribution functions for hadrons (/i), fq{x) stand for the distribution 
of the momentum fraction x carried by a parton inside a hadron. In analyzing the deep-inelastic electron-proton 
scattering, the parton distribution functions becomes crucial for instance. We note that these two nonperturbative 
functions are inter-related analytically via the Drell-Levy-Yan (DLY) relation [12]. Note that this analyticity is 
possible only if those functions can be described by the same function, defined in the different regions p!2] . 

Those functions have been studied intensively for several decades but still not fully-understood p!3H39] . Empirically, 
their information can be extracted from the available high-energy lepton-scattering data by global analyses with 
appropriate parameterizations satisfying certain constraints [13-18 . 

From theoretical points of view, there have been numerous works done for those functions so far: The momentum 
sum rules for the fragmentation functions were proven in terms of QCD in a rigorous manner [21 . In Refs. [23]426] , the 
Nambu-Jona-Lasinio (NJL) model was applied for these functions with the quark-jets and resonances, which satisfied 
the momentum sum rules. It turned out that the quark-jet contributions provide considerable contributions to the 
various fragmentation functions at the small z region. The Dyson- Schwinger (DS) method presented considerably 
successful results for the valance-quark distribution functions for the pion and kaon [27]. Monte-Carlo simulations 
with supersymmetric (SUSY) QCD were carried out to obtain the fragmentation function up to a very high energy in 
the center-of-mass frame ^/s The collins fragmentation functions which play an important role in the transverse- 
spin physics have also been studied in the quark-PS meson coupling model [35ti37]. Dihadron fragmentation functions 
were also studied in the same theoretical formalism [38] , 
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FIG. 1: Schematic figures for tfie quark-distribution (left) and fragmentation (right) functions, in which the soUd and dash 
hues denote the quark and pseudoscalar meson, respectively. 

Similarly, the quark-distribution function has also been studied with a nonlocal effective-chiral Lagrangian, em- 
ploying the DS method [28 . In Ref. [29 , the authors made use of the the Drell-Yan process including the soft-gluon 
resumption for the valance-quark distributions. A review for the experimental and theoretical status for valance-quark 
distribution for the nucleon and pion is given in Ref. [30 . 

In Refs. [31-34 , the instanton- motivated approaches were taken into account for computing the quark distribution 
functions, manifesting the nonlocal quark-pseudoscalar (PS) meson interactions. Employing the nonlocal chiral quark 
model (NLChQM) which is motivated by the dilute instanton-liquid model (LIM), which is properly defined in 
Euclidean space in principle j4Tti45] . NLChQM have accumulated successful results for various nonperturbative 
quantities in good agreement with experiments as well as lattice QCD (LQCD) simulations. [46H5T] . In [40] we use 
this model to calculate the fragmentation functions and the quark-distribution functions of the pions. There we limited 
ourselves in the SU(2) flavor sector. Our result is substantially different from the results of other models with local 
quark- meson couplings. In the present work, we extend our studies into the three flavors, {u^ s), i.e. those functions 
for the positively charged pion and kaon, i.e. 7r+ and taking into account the explicit flavor-SU(3)-symmetry 
breaking. We present the result with the two different explicit SU(3) flavor breaking patterns associated with two 
sets of the current quark masses. Furthermore, we figure out that our result of kaons fragmentation functions are 
significantly different from the ones from the other models. The calculated distribution and fragmentation functions 
are evolved to high by the DGLAP evolution using the code OCDNUM ^2\\53\ to compare with the empirical data. 
By doing so, we see qualitative agreement between our result and the empirical data, although some underestimates 
are shown in the small-(x,z) regions, according to the absence of the quark-jet contributions possibly [23-26 . 

The present report is organized as follows: In Section II, we briefly introduce NLChQM that we use and sketch our 
computation procedure. In Section III we present our numerical results and related discussions. The final Section is 
devoted for the conclusion and future perspectives. 



II. NONLOCAL PION-QUARK COUPLING 

The process that an off-shell quark (q) is fragmented into an unobserved set of particles (X) and on-shell hadron 
(/i), i.e. q ^ hX. is described by the unpolarized fragmentation function D^. A schematic figure for this function is 
given in the right of Fig. [l] It can be written with the light-cone gauge as follows [35] [36j : 

D'^iz, fe2 , ^) = j-J dk+Tv [A{k,p, m)7-] U-=p- , (1) 

where Note that indicates the flavor factor for the corresponding fragmentation processes. The values for for 
the PS meson, i.e. h = (p = (7r,i^), are given in Table [TT| in Appendix, /c, p, and z indicate the four- momenta of the 
initial quark and fragmented hadron, and the longitudinal momentum fraction possessed by the hadron, respectively. 
k± denotes {ko ± ks)/V2 in the light-cone coordinate. All the calculations are carried out in the frame where k± = 0. 
Here the 2:-axis is chosen to be the direction of k. On the other hand, kr = k — [{k ■ p)/|pp]p, defined as the 
transverse momentum of the initial quark with respect to the direction of the momentum of the produced hadron, is 
nonzero, /i denotes the renormalization scale at which the fragmentation process computed. Note that we consider 
this renormalization scale is almost the same with the momentum-transfer scale, i.e. ja^ ~ for simplicity, unless 
otherwise stated. The correlation A(/c,p, /i) reads generically: 

A(fc,p,M) = Y.J 0^e+"^<{om)\h,X){h,X\m\O)- (2) 
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Here represents the quark field, whereas ^ the spatial interval on the light cone. Furthermore one can integrate 
over fcr, 



D'1{z,ii) = ttz' 



Jo 



(3) 



The factor of z'^ is due to the fact that the integration is over = p — [[k - p)/\k\'^] the transverse momentum of 
the produced hadron with respect to the quark direction, and there is a relation between p± and /ct- P± = —zkT- 
The integrated fragmentation function satisfies the momentum sum rule: 



/ y2zD^{z,fi)dz = i, 

-Jo u 



(4) 



where h indicates for all the possible hadrons fragmented. Eq. Q means that all of the momentum of the initial quark q 
is transferred into the momenta of the fragmented hadrons. From the Drell-Levi-Yan (DLY) relation p!2 l [22l l23l [25l [36] , 
can be related to the parton distribution function provided that there is a proper analytic continuation. The 
relation is as follows, 



D'^^iz) = ^-f^ (x) , where x=-^, 



(5) 



where x denotes the momentum fraction possessed by a parton inside the hadron. A schematic figure for the quark- 
distribution function is depicted in the left of Fig. [l] 

In this article, we use NLChQM to investigate these nonperturbative objects, i.e. fragmentation and parton 
distribution functions. This model is motivated from the dilute instanton liquid model [4TH45]. We note that, to 
date, various nonperturbative QCD properties have been well studied in terms of the instanton vacuum configuration 
and the results are comparable with experiments and LQCD simulations |46ti49] . In that model, nonperturbative 
QCD effects are deciphered by the nontrivial quark-instanton interactions in the dilute instanton ensemble. However 
this model by nature is defined in Euclidean space because the (anti)instantons are well defined there by signaling 
the tunneling between the infinitely degenerate QCD vacua. Although there have been no rigorous derivation on the 
analytic continuation from the instanton physics to those in Minkowski one, there are still several challenging studies 
which try to apply the idea of the instanton physics to the physical quantities defined properly only in Minkowski 
space, such as the light-cone wave function [3TH3T. Following those studies we adopt the effective chiral action 
(EChA) from NLChQM in Minkowski space as follows: 



5eff [rriq 



-Sp In 



(6) 



where Sp and rhq denote the functional trace Tr J d'^x{x\ ■ ■ ■ \x) over all the relevant spin spaces and SU(3) current- 
quark mass matrix, diag(m^, m^^, mg), respectively. Note that, in deriving EChA in Eq. (|6|, we simply replace 
the Euclidean metric for the ( ant i) instanton effective chiral action into that for Minkowski space. The momentum- 
dependent effective quark mass generated from the interactions between the quarks and nonperturbative QCD vacuum, 
can be written in a simple n-pole type form factor as follows 



M{d^) = Mo 



nh? — d'^ -\- ie 



(7) 



where n indicates a positive integer number. We will choose n = 2 as in the instanton model [33l|34]. A stands for the 
model renormalization scale. It is related to the average (anti) instanton size p in principle. The nonlinear PS-meson 
field, i.e. U^^ takes a simple form [45 with the normalization following Ref. [35 to be consistent with the definition 
of the fragmentation function in Eq. |T]): 



/7^^((^) =exp 



^75 (A » 



1 



2Fs 



SFi 



where and A'^ stand for the weak-decay constant for the PS meson (j) and the Gell-Mann matrix. We note that, 
however, the value of F^^ can be determined rather phenomenologically to reproduce relevant physical quantities and 
conditions, even in NLChQM, and will discuss this in detail in Section III. We also write explicitly the flavor SU(3) 
octet PS-meson fields: 



%/6 



V 



K- 



K+ 



(9) 
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By expanding the nonlinear PS-meson field up to 0(^1) from EChA in Eq. one can derive the following effective 
interaction Lagrangian density in the coordinate space for the nonlocal quark-quark-PS meson vertex: 



— q^M{d^)^,{T-4>)^M{d^)q. 



(10) 



As we have done in our work in SU(2) sector [40 , we reach a concise expression for the elementary fragmentation 
function q (j)q' from NLChQM: 



MkMr 



M)) 



2MfM 



flVlf 



2k -p) 



8n^z{l-z) 2Fl 
here the momentum dependent effective quark mass is given by 

2A2 

Ml = Mo 



2A2 - £2 + 



(11) 



(12) 



and some notations are defined as follows: = ruf -\-M£ and Mf = ruf -\-Mo. Note that / and indicate the flavors 
for the initial {q) and final {q^) quarks, respectively. Note that the value of Mq can be fixed self-consistently within the 
instanton model [4TH49] with the phenomenological (anti)instanton parameters p ^ 1/3 fm and ^ ^ 1 fm, resulting 
in Mo ~ 350 MeV by the following self-consistent equation, defined in Euclidean {E) space. The masses for the pion 
and kaon are chosen to be m-j^^K = (140,495) MeV throughout the present work. Collecting all the ingredients above, 
one is led to a final expression for the elementary fragmentation function: 



[z^kl + [{z-l)Mj + Mf,]^] 



87r3 2F| [z^kl + z{z - l)M] + zM], + {I - z)mlf' 



(13) 



where and Mr are the momentum-dependent quark mass manifesting the nonlocal quark-PS meson interactions, 
read: 



[z'^kl + z{z - 1)(2A2 - (52) + zM} + (1 - z)mlY ' 



Mr 



Mo (2 A' 



2^2 



(2A2 -M2 )2- 



(14) 



As in [40] , a free and finite- valued parameter 5 in the denominator to avoid the unphysical singularities which appear 
in the vicinity of (z^kr) = 0, due the present parametrization of the effective quark mass as in Eq. At the 
renormalization scale in our model, the elementary fragmentation function is assumed to be able to be evaluated 
further by integrating Eq. (13) over kr as: 



(z, /i) = 27rz' 



f 

Jo 



(z, fi) kr dkT^ 



(15) 



where the factor z'^ in the right-hand-side is again comes from the integration over p± = —kr/z. Actually the 
connection between d^{z,p) and d^{z^k'^^ fi) would be far more complicated in principle [5^. We will present the 
numerical results for Eq. ([l5| in the next Section. 

Using the DLY relation in Eq. (jsj), one can derive the quark-distribution function as follows: 



fq{Xj M) 



3C,* MkMr 



[{x - l)Mf - xMf,f 



47r3 2F2 [kl + (1 - x)M'j + xMj, + xix - l)m2]2 ' 
where we have defined the effective quark masses for the quark distribution function by 

4MoA4(l - a;)2 , ^ 4MoA'^ 



Mk 



[fe2, + (1 - a;)(2A2) + xM], + x{x - l)m2]2 ' 



Mr 



(2A2 -M2^)2 



(16) 



(17) 



We see that Eq. ( |16| ) is equivalent to that given Ref. [24 , except for the momentum-dependent effective quark mass 
and qq<p coupling. Note that Al^ in Eq. (17) does not suffer from the unphysical singularities unlike that in Eq. (14), 



according to the different kinematic situations between those functions. Also it ensures the correct behavior that the 
quark-distribution function becomes zero as x ^ 1 [55^ . 
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FIG. 2: (Color online) Quark-distribution functions, multiplied by x, i.e. xff{x) for (0, Q') = (tt^u) (solid), (0, Q') = {K,u) 
(dot), and (0,^) = (^i s) (dash), at = = 0.36 GeV^. The thick and thin lines indicate the results from Model I and II, 
respectively. The curves for xfu{x) for Model I and II show negligible differences as shown. 



Similarly, the integration over kr can be performed as follows, resulting in a function of x at a certain renormalization 
scale /i ^ A: 



POO 

Jo 



(18) 



The minus-type quark-distribution function, which is nothing but the valance-quark distribution, satisfies the following 
normalization condition: 



/" 

Jo 



dx 



ftix,^i)-fiix,^i) 



(19) 



where Uq in the right-hand-side indicates the valance-quark number and becomes unity. We will tune the value 
to fulfill the normalization condition in Eq. (19) in the next Section. 



III. RESULT AND DISCUSSION 



In this section, we present and discuss our numerical results. We only compute the fragmentation functions for 7r+ 
and because the other PS-meson channels can be easily obtained by multiplying the flavor factors in Table \n\ It is 
due to the fact that we are only focusing on the elementary process here. For simplicity, we will denote the positively 
charged pion and kaon simply as tt and K hereafter. The values of F-j^^k in Eqs. (13) and (16) are determined by 
the normalization conditions for the quark distribution functions as in Eq. (19). In this phenomenological approach. 



we take into account two models with different patterns of the explicit flavor SU(3) breaking according to their 
current-quark masses: 



• Model I 

We consider the isospin symmetry for the u and d quarks, and 



Aqcd: {mu.md.rris) = (5,5,150) MeV. 



= 88.07 MeV, 



rpu- 



^Ks 



= 118.30 MeV, = 103.20 MeV. 



(20) 



• Model II 

We employ the current-quark masses from the PDG center values [56 : {mu^rrid^ms) = (2.5,5, 100) MeV. 



88.49 MeV, F^-^^"" 



116.60 MeV, 



95.94 MeV. 



(21) 
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FIG. 3: (Color online) Elementary (left) and renormalized (right) fragmentation functions, multiplied by i.e. 
zD'^(z), at — — 0.36 GeV^. The thick and thin lines indicate the results from Model I and II, respectively. 



and 



Note that F^^^^ and Ff^^^ are about 10% smaller than their empirical values, Ft^^k ^ (93, 113) MeV. This has 
already been observed in Refs. [40[|57] and explained by the fact that the so-called nonlocal contributions proportional 
to dM{p^) / d\p\ and the meson-loop contributions corresponding to large- A^c corrections have been neglected. Due 
to the similar reason, F'^~^^^ turns out to be slightly larger than its empirical value. Instead of including the 
nonlocal contributions in the present work, we have modified the F^^^k values as in Eqs. (20) and (21), satisfying 
the normalization condition for f^{x) for brevity [40]. This treatment is supposed to compensate the absence of the 
nonlocal contributions phenomenologically. Using the values in Eqs. ( [2Q| ) and ( [21] ), we first present the numerical 
results for the quark-distribution functions, multiplied by x in Fig.pl i.e. xf^{x) for (^,(7) = {tt^u) (solid), (^, ^) = 
{K^u) (dotted), and (^, ^) = (K^s) (dashed), at = = 0.36 GeV . The thick and thin lines represent the results 
from Model I and II, respectively. We find that the curves of xf^{x) from Model I and II are not distinguishable in 
the given resolution. It means the isospin-symmetry breaking effects are negligible as expected. Similarly the curves 
of xf^ {x) from different models show very small difference. In contrast, there is a clear difference between the two 
curves of xf^ {x) in Fig. In the other words, the fraction of momentum of the u quark is more sensitive to the s 
quark- mass than the one of the s quark inside the kaon. In terms of shape, we obtain very symmetric curves for xf^{x) 
around x = 1/2, whereas xf^g{x) are both tilted slightly to one side. For other applications such as the DGLAP 
evolution via QCDNUM17, which is a FORTRAN program that numerically evolves parton densities or fragmentation 
functions up to NNLO in perturbative QCD [52l|53], it is very useful to parameterize the quark-distribution functions 
in the form as 



xff{x) = A'l;x''t{l-xft^ 



(22) 



where A^, a^, and f3^ denote real constants. Using Eq. (22), we can parameterize the present numerical results shown 
in Fig. [2|as follows: 



Model I 



x/^(x) =3.22x^-^^(1 -x)^•^^ xf^{x) = 7mx^-^^{l-xf-^\ x/f(x) =8.77x^-^^(1 -x) 



1.74 



(23) 



• Model II 



xf^{x) = 3.21 x^-'^ (1 - xy-^^, xf^ix) = 8.51 x^-^^ (1 - xf-^\ xf^{x) = 9.97 x^-^^ (1 - x) 



1.84 



(24) 



After fixing F^ k for Model I and II, now we are in a position to present the fragmentation functions in Eq. (13). In 
the left panel of Fig. [3| we show the numerical results for the elementary fragmentation functions, multiplied by 
i.e. zd'^{z)^ at = A^ = 0.36 GeV^, in the same manner with Fig. [2J The thick and thin lines indicate the results 
from Model I and II, respectively. We observe that the pion fragmentation is much larger than those of the kaon. 
Furthermore, the curves of zd^ and zd^ are titled to the region z > 0.5. The differences between Model I and II for 
are all not significant in the case of fragmentation functions. 
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FIG. 4: (Color online) Plus-type (left) and min us-type (right) quark-distribution functions, multiplied by x, i.e. x[ff {x)^ff (x)] 
for = (7i^,K^) and q = {u,s), at ^ = 0.36 GeV^ (solid). The dot and dash lines indicate the LO and NLO DGLAP 
evolved functions, respectively, at — 4GeV^. The thick and thin lines indicate the results from Model I and II, respectively. 
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In Refs. [22^, 'SS^, it was argued that the elementary fragmentation functions do not hold the momentum sum rule in 
Eq. Q. To remedy this problem, the PS-meson cloud effects were devised and enumerated by the following relation: 

^ f\zdf{z)^M,, (25) 

where N'q stands for the renormalization constant representing the PS-meson cloud effects. The summation runs over 
all the possible PS-meson states which can be fragmented from q. As for the PS-mesons fragmented from one u 
quark, we have to sum over (I)' = (7r+, 7r°, K+) to obtain JVq. The values of Cf are listed in Table II in Appendix. Note 
that we do not include r] and rj' in cj)' because one needs calculate flavor-singlet contribution if one wants to include 
them in (j)^ and it is beyond our current framework. 

Typical value of A/" is estimated as (0.1 ^ 0.2) in other models. For instance, in according to the NJL model | 23l - 
[26] or quark-PS meson coupling models ptiSTj. it was given that AT^J^ = 0.03, A/"^^ = 0.08, and A/"^^ = 0.17, 
respectively. Here the superscripts PS and PV denote the pseudoscalar and pseudovector quark-PS meson couplings. 
Numerically, we have A/'n = 0.42(0.41) and ATg = 0.04(0.03) for Model I (II) in the present framework. Note that 

our value for JVg is underestimated because we ignore the contribution from d]'^ (z) here. These values are obviously 
larger than those in other models. It tells us that substantial PS-meson cloud effects have already been included in our 
model as discussed in the previous work [40 . Actually the momentum-dependent effective quark mass contains part 
of the PS-meson cloud effects, since in usual Dyson- Schwinger approaches [27 the effective quark mass corresponds 
partially to the dressed quark mass generated by the PS-meson could. Although the origin of the effective mass in 
our model is the nontrivial interactions between the quarks and (anti)instantons. Nevertheless, one still can argue 
that both of PS meson cloud and the quark- (ant i)instantons interaction are both from the QCD vacuum. Thus we 
define a renormalized fragmentation function for the pion and kaon as follows [40] : 

Dtiz) = (26) 



In the right panel of Fig. |3j we show the numerical results for zD^{z) via Eq. (26) in the same manner with the 
left panel of Fig. [sj The curve of zDf {z) is enhanced a lot compared to the curves of zD^{z) and zD^ [z) because 
Mu is much larger than Ms- zD^ is far more smaller than zD^ reflects the reality that the energtic u quark is most 
likely to be fragmented into pions rather than kaons. Moreover, we find the tendency zD^ <C zD^ for all the z 
regions. This large imparity is a quite different from the NJL- Jet model [24] where zD^ ~ zD^/2 approximately. 
Although they were computed at a different value from ours, = 0.2 GeV^ and their results have included 
quark-jet contributions which is supposed to be very important at low z regime. The difference of zD^ / zD^ between 
their results and ours at high z regime is still interesting. Presumably it is due to the nonlocality of the quark-meson 
interactions in our model. 

To compare our results with the empirical data which are usually extracted from high-energy scattering experi- 
ments [58], we need perform the high-Q^ evolution. To this end, we will employ QCDNUM17 [16, 52, 53 . As the 



inputs for QCDNUM17, we use the parameterized quark-distribution functions in Eqs. (23) and (24), taking into 
account the elementary processes only for the positively charged pion and kaon: 

xUZ - fl) = xfZ, xifS - /I) = -xfl xU: - m = 0, xfl = 0, x/l = x/5 = 0, 

x{f^ - fi) = xf^, x{f^ - /f ) = 0, x(/f - /f ) = -x/f , xf^ = 0, xff = 0, x/f = x/^, (27) 

where we have assumed that xf^ = x/J, due to the isospin asymmetry. In Fig. [4j we show the numerical results 

for the plus-type (left) and minus-type (right) quark-distribution functions, multiplied by x, i.e. x[f^{x) ±f^{x)] for 
(^, g) = (tt^u) (first row), {(j)^q) = {K^u) (second row), and {(j)^q) = {K^s) (third row). We note that the minus- type 
function is nothing but the valance-quark distribution function. The results at = 0.36 GeV^ are given in the solid 
line, whereas the dotted and dashed lines indicate the LO and NLO DGLAP evolved functions at = 4GeV^, 
respectively. The thick and thin lines indicate the results from Model I and II, respectively. The empirical data 
for the pion, given in the dashed lines, are extracted from a consistent NLO analysis of several high-statistics tt^N 
experiments including both Drell-Yan and prompt photon production [18]. 

As for the pion case. Model I and II provide almost the same curves as expected. As increases, the plus- 
type functions are tilted to the smaller x region due to the increasing gluon contributions. Although the plus-type 
function underestimates the empirical data in the vicinity of x = 0, the minus-type one reproduces the empirical data 
qualitatively well. The present underestimates in the small x region may be improved by including the quark-jet 
contributions as in the NJL- Jet model I24H26I. Here we do not consider those contributions. The second and third 



9 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.5 0.8 1 



z z 



FIG. 5: (Color online) The fragmentation functions zD^ (left) and zD^ (right) via LO and NLO DGLAP evolution at 
q2 ^ 1.0 GeV^ The empirical curves are taken from 14 . The thick and thin lines indicate the results from Model I and 11, 
respectively. 



rows show those functions for the kaon. We see again similar tendency with that observed in the pion case. The 
dependence becomes most visible for {(\)^q) = {K^u). 

In Fig. [sj we draw the renormalized fragmentation functions, multiplied by z and LO DGLAP evolved to = 
1 GeV^, for (^, q) = [K^ u) (left), and ((/), q) = [K^ s) (right). The dotted lines represent the empirical curves taken from 
a global analysis of charged-hadron {h) production data in electron-positron annihilation, i.e. e++e~ h^-X [14]. 
The shaded areas indicate the accumulation of the errors estimated by the Hessian method. In that analysis, the 
fragmentation functions for the PS mesons are parameterized as follows [14]: 

Ji/fcf) 

B(a5 + l,/3j + 2) 

where the NLO parameters, M^, a^, and are given in Table |l| B{a,b) in the denominator stands for the Beta 
function with the arguments a and b. The values for those NLO parameters are listed in Table [l] [14 . Comparing 
the empirical data with our results, we observe considerable overshoots around z = (0.4 ~ 0.5) exceeding the errors. 
Moreover, the results from Model I and II are quite close. 



In Fig. IIL we depict the minus- type quark-distribution functions multiplied by x, i.e x[f^ — f^], evolved to 
Q2 = 27GeV^ for the pion (left) and kaon (right). The LO (solid) and NLO (dotted) data for the pion is taken 
from Ref. [16] and are extracted from the muon-pair production experiment by 252 GeV pions on tungsten target. 
The Model I and II results are depicted by the thick and thin lines, respectively. As shown in the left panel, the 
isospin-symmetry breaking is negligible. It turns out that the numerical results reproduce the data qualitatively well. 
Note that similar results were already observed in Ref. [40 . In the right panel, we show the the kaon case fragmented 
from q = u for LO (solid) and NLO (dotted) and q = s for LO (dashed) and NLO (long-dashed) in the same manner 
with the left panel for Model I and II. The dependencies are obviously shown for the case for q = being similar 
to those in Fig. [4| The overall shapes of the curves for the kaon are very similar to those for the pion for the region 
X < 0.4, but the kaon curves drop more stiffly in comparison to the pion case when x > 0.5. 





Mt 


cA 






0.401 ±0.052 


-0.963 ±0.177 


1.370 ±0.144 




0.0740 ± 0.0268 


-0.630 ± 0.629 


1.310 ±0.772 


^ s 


0.0878 ± 0.0506 


2.000 ±2.913 


2.800 ± 1.313 



TABLE I: Parameters for Eq. (28) at = 1 GeV^, taken from a global x analysis of e+ ± e~ ^ /i ± AT [T] 
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FIG. 6: (Color online) Valance-quark distribution, x[ff(x) — f^{x)] for the pion (left) and kaon(riglit) via LO and NLO DGLAP 
evolutions at = 27GeV^. The experimental data of the ion is taken from 16 . The thick and thin lines indicate the results 
from Model I and II, respectively 



For simplicity we introduce a notation for the valance-quark distribution function as follows: 



(29) 



where q = {u^d^s) and (j) = {tt^K). In Fig. [7| we show the numerical results for uk/utt (left) and sk/u^^ (right) as 
functions of x. The results from Model I and II are depicted in the thick and thin lines, respectively. The results are 
evolved to = 27 GeV^ by the LO (solid) and NLO (dotted) DGLAP evolutions. The experimental data given in 
the left panel of Fig. [t] are taken from the 150 GeV incident-beam experiment for (i^~,7r~) + nucleus fi^ ii~ X j58] . 
A fitted curve is also given (dashed): ur/u^ = 1.1(1 -x)^-^^ ^ 

As for uk/u^, the numerical result reproduces the data qualitatively well. The curves decrease with respect to x 
when X > 0.4. The results of Model I show considerable underestimate for x > 0.5. As for the smaller in Model II, 
the numerical results are enhanced and become close to the data points. We also note that the LO and NLO results 
have negligible difference. It was argued that the value of uk/u^^ around x = 1 manifests the nonperturbative nature 
of the distribution function, since it is not affected by the high-Q^ evolution [27l [30l |55^ . In the literatures, the value is 
given as uk{x 1)/ut^{x ^ 1) = (0.2 ~ 0.5). From our numerical results for Model I, we have uk{^)/ut^{^) ^ 0.25, 
whereas it becomes for 0.31 Model II. Analytically, the quark-distribution function in the present theoretical framework 
at X = 1 reads from Eq. (I29|): 



f^{x ^ 1,/i) (X 



1 



2mqMo 



Mi 



(30) 



where q = (u^s) for (j) = (tt^K). Using the above equation and the values for Model I in Eqs. (20) and (21), we can 
make a rough estimation for the ratio as follows: 



UKjx ^ 1) _ F2 
Ut^{x 1) 



^ K 



1 



{rris - rriu) 
Mo 



0.32, 



(31) 



which is compatible with 0.25. Moreover, it is clear that the ratio in Eq. (31) is a function of typical nonperturbative 
quantities in QCD, such as the PS-meson weak-decay constants and constituent-quark mass Mq, in addition to 
the flavor SU(3) symmetry-breaking effects rriu <C m^. It is worth mentioning that different forms for the ratio were 
suggested as {F^/ Fk){Mu/Ms)^ from the BSE calculation [27 and {Mu/MsY from the NJL calculation [55 . In the 
right panel of Fig.[7| the ratio sk / '^tt is given as a function of x. The behavior of the curves are very different from those 
for uk/ut^ which can be understood easily from Fig. Ill The value at x = 1 becomes 1.18 and 1.01 for Model I and II, 
respectively, which is about (3 ~ 5) times larger than that for uk/ut^- This larger ratio is caused by that the second 
term in the square bracket in Eq. (31 ) goes to zero for s^/'^tt, i-e. q = uioi the both of the valance-quark distribution 
functions. Approximately, using Eqs. (31), (20), and (21) for Model I, we obtain sk{x 1)/ut^{x ^ 1) ~ 0.84, which 



is still compatible with the numerical result. 
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FIG. 7: (Color online) Ratios of the valance-quark distributions, uk/uj^ (left) and sk/uj^ (right) via LO and NLO DGLAP 
evolutions at = 27GeV^. In the right panel, we show the experimental data from Ref. [58J and the fitted curve (dash) from 
Ref. 23- The thick and thin lines indicate the results from Model I and II, respectively. 

IV. SUMMARY AND OUTLOOK 

We have investigated the fragmentation functions of the pion and the kaon employing NLChQM in the light-cone 
coordinate and apply the DLY relation to extract the quark-distribution functions from the fragmentation functions. 
Two patterns of SU(3) flavor breaking have been considered: one is for the current-quark mass {mu^md^rris) = 
(5, 5, 150) MeV (Model I) and another is for (m^, m^, m^) = (2.5, 5, 100) MeV (Model II). In what follows, we list the 
important procedures and observations in the present work: 

• Ft^.Ki are determined by the normalization conditions for the valence-quark distribution functions. By doing so, 
we compensate the absence of the nonlocal contributions for the vector-current conservation phenomenologically. 
The calculated values of Ft^^k turn out to be about (5 ~ 10)% smaller from their empirical values. 

• The curves of the quark-distribution function, multiplied by x, i.e. xff{x) for the pion and kaon exhibit very 
different shapes. It is due to the explicit flavor SU(3) symmetry breaking. Moreover, the shapes of those curves 
are substantially different from the ones generated from usual local-interaction models. 

• The curves of the elementary fragmentation functions of the pion and kaon, multiplied by i.e. zd^ are also 
different. The pion one is symmetric around z = 1/2 but the kaon ones are titled to higher z side. Furthermore 
the magnitude of the pion one is large more larger than the kaon ones. 

• For the renormalized fragmentation functions of the pion and kaon, multiplied by i.e. zD^, we find zD^ are 
much larger than zD^ . This is quite different from the results from other models with the local couplings. 

• The kaon fragmentation functions are evolved to = 1 GeV^. Compared with the empirical curves our curves 
are much larger. We also notice that zD^ is more sensitive to the strange quark mass than zDf . 

• The plus- type and minus- type quark-distribution functions are evolved to = 4 GeV^. They are compared with 
the empirical data for the pion, resulting in a qualitative good agreement with them. We also computed those 
for the kaon within the same framework. The fragmentation functions for the kaon are evolved to = 1 GeV^ 
and we observed that strong overshoots in the vicinity of 2; = 0.5 compared with the empirical data. 

• We presented the valance-quark distribution functions at = 27GeV^ for the pion and kaon. The numerical 
results are qualitatively compatible with the empirical data for the pion. We also calculated the ratios for the 
valance-quark distributions for the kaon and pion. Although the numerical results showed qualitative agreement 
with the data, there appears considerable deviations for the region x > 0.5 for uk/u^- Theoretical estimations 
for sk/utv are presented. 
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As noticed in Section III, we have not taken into account the axial-current conservation in the present framework, 
which may become problematic for the nonlocal quark-PS meson interactions as in NLChQM [33 . Moreover, we also 
note that, in the present framework, the momentum-dependencies in the effective quark masses in the fragmentation 
and distribution functions, given in Eqs. (13) and (16), were simplified to avoid complexities in numerical calculations, 
i.e. Ml Mq. In order for theoretical consistencies, these two issues might be taken into account carefully, although 
the present results are phenomenologically acceptable by comparison with the data. We verified that the (valance) 
quark-distribution functions are considerably modified by correcting those issues, satisfying relevant theoretical con- 
straints, and the same for the fragmentation functions accordingly [59 . Related works are under progress and appear 
elsewhere. 
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TABLE II: Flavor factors in Eq. ([T}. 
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